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Abstract 

Magnetic phase transition (MPT) to magnetic quasi-long-range order (QLRO) phase in a three- 
dimensional Heisenberg weak (D/J = 4) random anisotropy (RA) model is investigated by Monte 
Carlo simulation. The isotropic and cubic distributions of RA axes are considered for simple- 
cubic-lattice systems. Finite-size scaling analysis shows that the critical couplings for the former 
and latter are K^. = 0.70435(2) and K^. = 0.70998(4), respectively. While the critical exponent 
1/z/ = 1.40824(0) is the same for both cases. A second-order MPT to the QLRO phase is therefore 
evidenced to be possible in favor with the existence of the QLRO predicted by recent functional 
renormalization group theories. 
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Rare-earth amorphous magnets are among magnetic glasses where the randomly- 
quenched disorder is in the form of intrasite random anisotropy This characteristic 
was first realized by Harris, Plischke, and Zuckermann in the Hamiltonian of their random 
anisotropy model (RAM) 2|, 

H = - J ^ 5, ■ 5, - D ^ (a, ■ Sif, (1) 

where Si is a unit vector representing the Heisenberg spin on site i of a simple cubic lattice 
of size L, which is subject to a local onsite anisotropy of strength D > and uniaxially-easy 
axis represented by a randomly-quenched unit vector Oj, and which can interacts with its 
six nearest neighbors 5*^ via ferromagnetic exchange interaction of strength J > 0. The 
degree of such a topological disorder is characterized by the disorder strength D/ J and the 
distribution of the easy axes p(aj) which is either of an isotropic form 

1 r(3/2) . , 

where T[x) is Euler gamma-function, or of an anisotropic (cubic) form 
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where S{y) is Dirac (5-function and are Cartesian unit vectors. Vector di is equally likely 
to point in any direction of the space for the former case while it points with the same 
probability of the latter in one of six directions along three Cartesian coordinate axes. 

RAM is the prototypical model system of a majority of amorphous magnets whose mag- 
netic properties have been intensively investigated since around 1973. The most interesting 
questions about the model to which answers are still not completely clear or are contra- 
dictory concern (i) the nature of magnetic phase ordering at low temperatures and (ii) the 
nature of magnetic phase transitions (MPTs). According to Imry-Ma argument 3|], long- 
range order (LRO) is destroyed in systems with a continuous symmetry in d < 4 spatial 
dimensions by even arbitrarily weak disorder of the form of either random field (RF) 
or random anisotropy (RA) [4']. Amorphous magnetic ordering (AMO) may be intuitively 
described by one of the non-collinear spin structures (NCSSs) suggested by Coey j^. The 
feature of the NCSSs is the possibility that spins within a domain may be frozen into more 
or less random orientations. Chudnovsky et al. [6,] proposed a phenomenological theory 

1 



based on the Hamiltonian in Eq. ([T]) for the description of AMO. Their theory shows that 
the characteristic of AMO depends crucially on the parameter A = {D/J){Rc/a)'^, where 
Rc is the scale of the spatial correlation of the easy axes and a is the atomic spacing. Of 
special interest is the prediction of the correlated spin glass (CSG) phase for weak RA, i.e. 
A <^ 1. The CSG exhibits a smooth rotation of the magnetization over the volume so that 
the directions of the magnetization are ferromagnetically correlated on a quite large length 
Rb > -Rc, (e.g., -Rd ~ A-'^Rc = {J/Dfa^/Rl for d = 3). Ahhough the ferromagnetic 
ordering exists on the scale i?D there are no sharp boundaries between the ferromagnetic 
regions, i.e. no sharp domain walls. In zero field, the net magnetization of the CSG is 
zero and the susceptibility is very large but finite x = |(^/4)~^ ~ {J / DY{a/ R^Y . For 
strong RA, when A ^ 1, a spin-glass-like (SG) state, the sperromagnet (SPM), is found 
with Rd ~ i?c- Recent theoretical work by Feldman [3] using the functional renormalization 
group (PRC) method has shown that in many cases LRO is prohibited in such the systems 
of continuous symmetry as RF and RA glasses, but instead quasi-long range order (QLRO) 
can emerge in the d = A — e dimensions for the presence of weak disorder. The QLRO is 
peculiarly characterized by a power-law correlation function and the average value of the 
order parameter (e.g. the net magnetization) over the volume being zero in zero field. The 
QLRO is more common in impure systems and not prohibited by non-Abelian symmetry. It 
is worth noting that there seems to appear a disagreement between the CSG and the QLRO 
theories that the CSG indicates large but finite correlation length and susceptibility while 
they are infinite for the QLRO. In our opinion, this paradox may be removed by the fact 
independently predicted in the theories by Cochrane et al. js] and by Elsasser et al. 9] 
that no correlation indeed exists between the easy axes of neighboring sites. One, therefore, 
should take the limit Rc ^ for Rd and x in the CSG theory above to approach the infinity. 

QLRO was clearly evidenced in a number of numerical experiments using optimal Monte 
Carlo (MC) methods 10, 11|. In particular, Itakura [ll^ has recently investigated the 



3D Heisenberg RAM in Eq. ([T]) with the easy-axis distribution in Eq. ([2]) and found 
that the spin-spin correlation function for the low temperature phase can be described by 
G{r) oc r~''~^ exp(— r/,^) with a nonuniversal exponent rj. The inverse of the correlation 
length ^ is finite for large values of D/J and vanishes when D/J< 5 so that a nonuniversal 
QLRO emerges with the correlation function of the frozen power-law form G{r) oc r~^~^. 
While a second-order MPT to the QLRO for a weak RA case D/J = 1 of the 3D XY RAM 
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ler the crossover from the paramagnetic 



ll| for a weak RA case D/ J = 4 of the 



was found by RoBler using MC simulations [12] , whet 
state to the QLRO, which was observed by Itakura 
3D Heisenberg RAM, is a true second-order MPT is still an open question. 

In this paper, we present our results of the MPT and critical behaviors for the case 
D/ J = Aoi the 3D Heisenberg RAM in Eq. ([T]) similar to that of Itakura's work. Systems of 
the simple cubic lattice with the periodic boundary condition are investigated. We, however, 
consider both cases of the distribution of the easy axes as in Eqs. ([2]) and ([3]) which are 
called isotropic radom anisotropy model (IRAM) and anisotropic radom anisotropy model 
(ARAM), respectively. The reason for this consideration of the easy-axis distributions shall 
be cleared up in the sequel. Since we do not study the glassy phase transition but the MPT 
we, therefore, focus only on the magnetic order parameter, that is the net magnetization 
m = ^MJTMJTMF where = ^. Si^. 

We first aim to determine the critical temperature, T^, and the critical exponent of the 
correlation length, u. We perform our important-sampling MC simulations with the single 
histogram method (SHM)[l^, [l^. Technically, it has been already optimal to use cluster- 
flip MC simulations (CFMC) and the multiple histogram method (MHM) to determine 
Tc and u independently in the pure Heisenberg model {13! . [l^ . Unfortunately, neither 
CFMC nor MHM can be applied efficiently for systems of the RA quenched disorder. The 
reason for the CFMC case has been pointed out by RoBler [l2|. For a given couple of 
lattice sizes, the intersection of their corresponding Binder's cumulant curves varies strongly 
when changing the set of quenched disorder variables {ai}^^^. Consequently, it is expensive 
to carry out several runs with various sets of {aj}^^ at different temperatures. Instead, 
we use the finite size effect (FSE) analysis for the maximum temperature, Tc{L), of the 
susceptibility x for each lattice size in the range 18 < L < 40 to determine simultaneously 
Tc and u with moderate efforts. Our simulation detail is as follows. For each lattice size, 
we run simulations for a modest number of equilibrium configurations (ECs) at different 
temperatures with the same fixed set of {aj}^]^ to calculate x(T) as a function of the 
temperature in order to roughly estimate Tc{L). We then carry out only a single long-run 
simulation at the estimated Tc{L) and store 10^ ECs 10 MC sweeps (MCS) apart. These 



stored ECs are subjected to the single histogram sampling as in [13|, [ij] to calculate x{T) 
and then to determine precisely Tc{L) for each set of {ctj}^!. Because of the nature of the 
SHM, we prefer the /T- representation for the sake of convenience, where K = J/T is called 
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the temperature coupling. Fig. [T] presents the size dependence of Kc{L). Each data point 
in the figure is averaged over ten independent sets of {ai}fli. Apparently, the correction-to- 
finite-size formula Kc{L) = Kc + aL~^^'^ describes quite well the asymptotic behavior of FSE 
without the need of the additional correction term L~'^ for both IRAM and ARAM cases. 
The best fit of this formula to the calculated data yields the values of the critical coupling 
(temperature) = 0.70435(2) (T^ = 1.41974(4)) and = 0.70998(4) (Tc = 1.40849(0)) 
for the former and the latter cases, respectively. While the value of the critical exponent of 
the correlation length is the same, i.e 1/z/ = 1.40824(0) or u = 0.70998(0), for both cases. 

Furthermore, we calculate Binder's cumulant Ul{K) = 1 — {m'^{K)) /3{m'^{K)y as a func- 
tion of K for lattice sizes 4 < L < 24 from the simulations of lO'' statistically independent 
ECs at various temperature couplings about K^. The calculated data, which are averaged 
over a large number of sets {aj}^^, and their rescaled data U[L^^'^{K/Kc — 1)] using the 
corresponding values of and 1/z/ above are shown in Fig. [2l It is remarked from the 
figure that for either case of the easy-axis distributions there is a well-defined intersection 
of the family of Binder's cumulant curves which indicates the fixed point of a second-order 
MPT. The Binder's cumulant, as a thermodynamic function, obeys the same finite-size scal- 
ing theory developed by Landau and Binder [13] as it does in the pure Heisenberg model. 
Strikingly, the rescaled Binder's cumulant curves for both IRAM and ARAM collapse nicely 
into a universal curve shown in Fig. [31 This fact means that systems of both types of 
the distributions belong to the same universality class with the correlation length's critical 
exponent u = 0.70998(0). 

n 

Since the work by Aharony pj| about the 'runaway' behavior, i.e. no second-order 
transition in the RAM in Eq. ([1]) for any dimension d < 4, the nature of phase transition of 
the model has become a longstanding question. Holovatch et al. have more recently showed 
theoretically that the second-order MPT is destroyed by the RA with the distribution in Eq. 
([2]) but survives in RAM with the distribution in Eq. ([3]). However, Monte Carlo simulations 
by RoBler 12| , by Itakura HI] and in our present work apparently contrast to the conclusion. 
In particular, we find nothing different in terms of the critical behavior between IRAM and 
ARAM as shown above. It is concluded that our present study together with RoBler's work 



121 ] shed the light to the possibility of a second-order MPT in the 3D RAMs of weak RA, 
which is in favor with the existence of the QLRO in 3D weak RAM predicted in the recent 
FRG theory by Feldman Q]. We hope that our work will stimulate further Monte Carlo 
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simulations on issues of the QLRO and the phase transition to it. 
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FIG. 1: Plots of KciL) vs. L-^/" for IRAM (circle) and ARAM (square). 
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FIG. 2: Binder's cumulant curves Ul{K) for IRAM (a) and ARAM (b). The insets show that the 
rescaled Binder's cumulant curves overlap into a universal curve for either case. 
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FIG. 3: The collapse of rescaled Binder's cumulant curves for both IRAM and ARAM. Systems of 
IRAM and ARM for the case D/ J = A may belong to the same universality. 
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